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Introduction 


As an analytic solution to the heavy-ion transport 
equation in terms of Green’s function representing 
nuclear and atomic/molecular processes, a heavy-ion 
transport code including a database has been pro- 
vided for laboratory ion beam applications. Results 
based on the new code were compared with pertur- 
bation theory results (ref. 1), which previously had 
been compared with those of 20 Ne transport exper- 
iments at the Lawrence Berkeley Laboratory (LBL) 
BEVALAC accelerator. (See refs. 2 and 3.) In the 
LBL comparison, the primary errors in the computa- 
tion were attributed to the nuclear cross sections and 
the approximations used in applying the acceptance 
functions. (See ref. 1.) The perturbation code was 
converted to access the NUCFRAG database (refs. 4 
and 5), then a direct comparison between the pertur- 
bation code and the nonperturbative Green’s func- 
tion code was made. (See ref. 1.) In this compari- 
son, the sequence of perturbation terms appears to 
be converging toward the nonperturbative result even 
though the nonconvergence of the lighter fragments 
in the first three perturbation terms is clearly ap- 
parent. (See ref. 1.) Although the nonperturbative 
Green’s function code eliminates the need to control 
truncation and discretization errors, it requires fur- 
ther development for comparison with space radia- 
tion codes. (See ref. 6.) 

Aside from the questions of numerical conver- 
gence and convergence of the series solution, the solu- 
tions themselves must represent the fields associated 
with all the isotopes produced in the fragmentation 
process. In principle, several hundred such isotopes 
would be required for the transport of iron beams. 
In practice, only a hundred or so such isotopes con- 
tribute to the solution in a significant way. The 
cosmic ray code that models high-charge and high- 
energy ions (HZE) and transport (HZETRN) devel- 
oped at Langley Research Center uses 59 isotopes, 
whereas many other simulations use only 29. (See 
refs. 6-8.) Determination of the number of isotopes 
required for an adequate laboratory beam simulation 
is the purpose of the present study. In this report, we 
recall the solution procedures and examine the effects 
on solution accuracy of representing a reduced set of 
isotopes. Although we consider only iron beams on 
epoxy targets, the conclusions are presumed applica- 
ble to other target materials and ion beams lighter 
than iron. This interpretation can be made because 
the isotopes produced are characterized by the nu- 
clear properties of the fragmented beam particles. 


Green’s Function for a Single Medium 

We restrict our attention to the multiple charged 
ions for which the Boltzmann equation may be re- 
duced (ref. 6) to 

L~WE^ {E) + aj ^( x - E ) = Y. a J k<Pk{x - E '> (!) 

-I k 

where <pj (x, E) is the ion flux at x with energy E 
(MeV/amu), Sj(E) is the change in E per unit 
distance, oy is the total macroscopic reaction cross 
section, and a y*. is the macroscopic cross section for 
the collision of ion type k to produce an ion of type j . 
The solution to equation (1) is found subject to the 
boundary condition 

<f> j (0,E) = f j (E) (2) 

For this boundary condition, laboratory beams have 
only one value of j for which fj{E) is not zero, fj (E) 
is described by a mean energy E 0 , and the energy 
spread a is such that 

fj(E) = -2=- exp [-(E- E 0 f/2a 2 } (3) 

v Z7TCT 

The usual method of solution is to solve equation (1) 
as a perturbation series. (See refs. 1 and 6.) In 
practice, the computational requirements limit the 
usefulness of the technique for deep penetration. (See 
ref. 3.) 

Green’s function is introduced as a solution of 

tc~W: ' S i {E . )+ °i C]m{x - E - Eo) = E E - Eo) (4) 

J k 

subject to the boundary condition 

G jm (0,E,E o ) = 6 jm 6(E-E o ) ( 5 ) 

The solution to equation (1) is given by superposition 
as 

<j> 3 {x y E) = [ G jk {x,E,E f )f k {E')dE' (6) 

k J 

If Gjh(x, E , E f ) is known as an algebraic quantity, the 
evaluation of equation (6) may be accomplished by 
simple integration techniques and the associated er- 
rors in solving equation (1) numerically are avoided. 
(See ref. 9.) 



The above equations can be simplified by trans- 
forming the energy into the residual range as 


r ,-fJ 

1 J Sj(E') 
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and defining new field functions as 

rpj(x,rj) = Sj(E)<t>j(x,E) 
G jm (x, Vj , r' m ) = Sj{E)Gj m (x, E, E' 
fj( r j) = Sj(E)fj(E) 
Equation (4) becomes 
d d 


upper and lower limits of equation (17) are switched. 
The higher terms are approximated as 

53 im 
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In the above expression 
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(9) 

(10) 
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and 
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with the boundary condition 

Gjm(0’ r j>r' m ) = Sj m 6(rj - r m ) (12) 

and with the solution to the ion fields given by 


(i) 

dent on x for i > 0, which we represent as Qj 7n (x). 
(See ref. 3.) In terms of the above arguments, the 
solution to equation (1) becomes (ref. 3) 


>rj) = e ajX fj( r j +«) 




( 21 ) 
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Note that i/j is the range scale factor as vjrj = v m r m 
and is taken as Vj — Zj /Aj. The solution to equa- 
tion (11) is written as a perturbation series 

Gj m (x, rj,r' m ) = 53 Gjll(x, rj,r' m ) (14) 


In equation (21), r f ml and r f mu are given by the lower 
and upper limits of the inequality in equation (17). 
The symbol F m (r f m ) refers to the integral spectrum 


lx; 
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We note that 


where 

= 9(])bjmt(x + Vj - r' m ) 

and 


F m {r' m ) = F m {E') 
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and 
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where ^^(x, rj, r f m ) is zero unless 


r ' m Is 
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V; , , Vj 

T -( r J + x ) < r m < ~r r J + x 


(17) 


for v m > v r If Vj > v m , as can happen in neutron 
removal, the negative of equation (16) is used and the 


We now introduce nonperturbative terms for the 
summation in equation (21). 

First, we recall that the g function of n arguments 
was generated by the perturbation solution of the 
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transport equation after neglecting ionization energy 
loss (ref. 1), which is given by 

d \ 

7^ = ^ ^ &jk9km{ x ) (26) 


subject to the boundary condition 


9jm (6) &jm 

(27) 

for which the solution is 


9jrn(x) = 6j m g(m) + cr jm g(j,m) + • • ■ 

(28) 

The equation 


9jrn(x) = J^9jk( x ~ y)9km(y) 

(29) 


k 


is also true for any positive values of x and y. 
Equation (29) may be used to propagate the function 
gjm{x) over the solution space, after which 

Gjm(x,rj,r' m ) « e~ a i x 6 jm 6{x + rj - r' m ) 

+ Vj [ffjrn(^) ~ e~ <TjX ^jm] ^ 
x{v m - Vj) 

The approximate solution of equation (1) is then 
given by 

Vj{x,r :) ) - e~ (T J x f(r J + x) 

E ^j \jHjm (^) “ C J 

x(v m - Vj) 
m 

x [E m (r mw ) — F m (r ml )J (31) 

Note that the computational procedures are affected 
by the size of the number of the elements in equa- 
tions (28) and (29). The number of terms in the 
application of equation (29) increases as TV 2 , where 
N is the number of isotope fields represented in 
the solution given by equation (31). For computa- 
tional efficiency, the goal is to minimize the number 
of isotopes without greatly compromising solution 
accuracy. 

Comparison of Isotope Tables 

The nonperturbative method generates Green’s 
function for any ion of charge Z < 28 that results 
from the impact of that ion on a material medium, 
including the secondary fragment fields. The atomic 
weight (taken as the nearest integral value) and 



Figure 1. Convergence of charge distribution of 505 MeV /amu 
iron beams in 5 g/crn 2 of epoxy for 59 and 80 isotopes. 

charge associated with each field function are trun- 
cated to the nearest isotope .4/, Z[ in the isotope 
table. The truncation minimizes the distance to the 
nearest isotope using the square-distance function 

D il = (A i -A l f + 4(Z i -Z l ) 2 (32) 

where A u Z } is the isotope produced in the fragmen- 
tation event. Clearly, the accuracy in the transport 
result requires the isotope list to contain the main 
isotopes produced in the fragmentation event; the 
isotopes of lesser importance may be approximated. 
Initially, 59 isotopes were selected to represent each 
nuclear mass value between 1 and 58. Such a list was 
found adequate for the transport of galactic cosmic 
rays using the HZETRN code. (See refs. 5 and 10.) 
However, such a representation was inadequate for 
the transport of an iron beam using the non- 
perturbative code GRNTRN (ref. 1); thus, more iso- 
topes were added to the table. The total flux of 
identified projectile fragment nuclei between H and 
Fe was found for 505 MeV/amu monoenergetic 56 Fe 
beams incident on epoxy of 5 g/cm 2 . The resin se- 
lected is tetraglycidyl 4,4 / diaminodiphenylmethane 
(TG 4, 4 f DDM) epoxy cured with diaminodiphenyl- 
sulfone (DDS). A repeat cured unit of epoxy molec- 
ular structures is C37H42N4O5S and has the density 
1.32 g/cm 3 . Epoxy was applied because it is a more 
common material and may be fabricated and sup- 
plied as a shield medium. The results based on tables 
of 59 and 80 isotopes are shown in figures 1 and 2. 
The integral output spectra for the projectiles and 
fragments in figure 1 show a somewhat similar charge 
distribution for both 59 and 80 isotopes, but the mass 
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ire 2. Convergence of mass distribution of 505 MeV/amu 
iron beams in 5 g/crn 2 of epoxy for 59 and 80 isotopes. 


Isotopes 



Figure 3. Convergence of mass distribution of 505 MeV/amu 
iron beams in 5 g/cm 2 of epoxy for 59, 80, 85, and 
100 isotopes. 

distributions in figure 2 show large differences. The 
80-isotope table is probably adequate for applications 
in which charge is the dominating factor (e.g., linear 
energy transfer), but the mass distribution could be 
substantially improved through an expanded isotope 
list. 

A modest change to the 80-isotope list was made 
with the addition of 5 isotopes; the results of this 
change are shown in figure 3. Significant improve- 



Figure 4. Convergence of charge distribution of 505 MeV/amu 
iron beams in 5 g/cm 2 of epoxy for 59, 80, 85, and 
100 isotopes. 



Figure 5. Convergence of mass distribution of 505 MeV/amu 
iron beams in 5 g/cm 2 of epoxy for 100, 108, and 
110 isotopes. 

merit in the mass distribution is achieved for A/ < 40, 
but the result degrades at higher mass numbers. The 
charge distribution was less accurate arid the use of a 
100-isotope list could not adequately resolve the con- 
vergence problem for the mass distribution as seen 
in figures 3 and 4. The isotope tables were incre- 
mentally expanded with continuous improvement in 
the mass distribution, as seen in figures 5 8. The 
final list of 122 isotopes appears to be the mini- 
mum set required to represent the fragment mass 
distribution. The charge distribution had nearly con- 
verged at 80 isotopes and no substantial change in its 
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igure 6. Convergence of mass distribution of 505 MeV/amu 
iron beams in 5 g/cm 2 of epoxy for 110, 113, and 
116 isotopes. 



Figure 7. Convergence of mass distribution of 505 MeV/amu 

iron beams in 5 g/cm 2 of epoxy for 119 and 122 isotopes. 

convergence occurs beyond 100 isotopes, as seen in 
figure 9. 

Although the specific tests were derived for an 
iron beam on a given epoxy resin, the isotope distri- 
butions are largely dominated by the nuclear physics 
of the projectile fragments and virtually all elements 
are produced below the projectile atomic number; 
thus, we expect similar convergence properties for 
other shield materials. We also note that the iron 
beam is a principal contributor to galactic cosmic 


Figure 8. Convergence of mass distribution of 505 MeV/amu 
iron beams in 5 g/cin 2 of epoxy for 122 and 125 isotopes. 
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Figure 9. Convergence of charge distribution of 505 MeV/ 
amu iron beams in 5 g/cm 2 of epoxy for 100 isotopes and 
greater. 

ray exposure and the current results indicate that 
the space shield calculations require a larger table 
than the 59 isotopes currently listed. 

Discussion of Results 

Computational precision of charge and mass dis- 
tribution is provided by adding more isotopes to 
the table. However, each additional isotope requires 
additional computation time to generate the non- 
perturbative Green’s functions for the selected iso- 
tope table. An optimal choice of isotope table is 
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Figure 10. Mass distribution of 505 MeV/amu iron beams in 
5 g/cm 2 of epoxy for revised table with 122 isotopes. 


needed for both computational precision and prac- 
tical computation time. 

The detailed isotope selections are shown in ta- 
bles I and II. Figures 8, 9, and 10 show that an opti- 
mal choice is the revised table of 122 isotopes. (See 
table II and attendant footnote.) In figures 8 and 10, 
the integral output spectra converge within 5 percent 
compared with the 125-isotope list over the entire 
projectile fragment nuclear mass range of 1 to 56. 
In figure 9, integral output spectra are plotted over 
the charge range of 1 to 26, where the maximum 
difference is within 2.7 percent over the entire pro- 
jectile fragment nuclei range for tables with 100 or 
more isotopes. A similar plot for 80 isotopes (fig. 1) 
gives a maximum difference of 3.1 percent. Although 
the largest list considered (125 isotopes) may not be 
fully converged, we believe that the error introduced 
by the 125-isotope list is much less than 5 percent in 
mass and 2 percent in charge distribution. 

Concluding Remarks 

Improvements in the treatment of the nuclear 
database are required so that space radiation codes 
will agree well with experiments. The improvement 
addressed in this research was the determination of 
an optimal isotope table to generate the nuclear data- 
base that gives both computational precision and 
practical computation time. An iron beam in epoxy 
was chosen to study the effects of isotope list selection 
on the mass and charge distributions of the trans- 
mitted fluence computed by nonperturbative meth- 
ods in the transport of high-charge and high-energy 
ions. A table of 80 isotopes gives charge versus flu- 


ence spectra that converge within 3.1 percent; a table 
of 100 isotopes converges within 2.7 percent. A ta- 
ble of 122 isotopes gives nuclear mass versus fluence 
spectra that converge within 5 percent. These tables 
also result in practical computation times. Iron is the 
most abundant massive ion in space and the fragmen- 
tation event is dominated by the nuclear structure of 
the projectiles, so these results are generally appli- 
cable to other materials and ions important to the 
space radiation problem. 
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